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Abstract 


Efficient control laws that ensure fast and accurate maneuver of Spacecrafts 
along with accurate and reliable control of the attitude are crucial for Space 
applications. In addition, the angular path traversed in such a maneuver is 
also an important aspect as this has a direct effect on the time as well as en- 
ergy required during such maneuvers. In this thesis, the structure of Control 
laws necessary for a shor test-path Eigen-axis rotation between two attitudes 
arc investigated. The attitude feedback taken is in terms of quaternion pa- 
rameters which is a represention of attitude in 3-Dimentional space well sujted 
for space-applications. Based on this, few Control laws are suggested. Global 
asymptotic stablity of the resulting close-loop systems are proved through Lya- 
punov stability analysis. Computer simulations of the proposed Control laws 
arc carried oiit to depict the responses of the syst.cms under various initial 
conditions. 
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Chapter 1 


Introduction 


In- space-technology applications, the attitude of a spacecraft in orbit is often very 
important. For example, communication satellites may have to be pointed in a par- 
ticular direction for recieving and beaming signals efficiently. Imaging satellites may 
be required to point accurately at some area on earth. For applications like space- 
telescopes, accurate steady pointing towards an object in space, for long periods of 
time, is absolutely crucial for observation and long-exposure photography. In appli- 
cations like space-based defence, fast and accurate pointing and tracking of targets 
become extremely critical. In addition, in recent years, there are also increasing 
applications of giant robot arms in space for repairing, re-orienting, and retreiving 
satellites. 

Evidently, it becomes extremely important to have a good control over the atti- 
tude of the spacecraft. In addition, other desired criteria like low fuel consumption, 
fast response, etc. make the attitude control of spacecrafts in orbit an important and 
challenging area of control applications. 

Attitude control of spacecrafts may be broadly categorised into two types. In some 
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cases, the body is desired to maintain a particular orientation in space. The control 
problem thus consists of bringing the body from any arbitrary initial orientation to 
this particular orientation. This class of control systems are termed as regulators, 
or of the regulatory type. On the other hand, in many applications the spacecraft 
is desired to move in a particular manner tracing a pre-defined trajectory. Here the 
control problem involves taking the body along this trajectory. This class of control 
systems are termed as those of the tracking type. 

In addition, another area in which the orientation of solid body in 3-Dimentional 
space is sought to be controlled, is robotics. In robotics, robotic hands are used for 
working tools, as in assembly lines, in handling materials, and similar cases. In these 
areas, often the fingers of the robot-arm act like actuators in controlling attitude of 
another rigid body, like a tool or vessel. Due to the inherent similarity of the control 
problems involved, many techniqes, and developments in the control of spacecrafts 
can be applied to robotics. 


1.1 Representation of attitude and its control 

For the purpose of analysis of general rigid-body motion, two reference-frames are 
considered, each consists of a set of three right-handed cartesian axes. One of them 
is fixed on the rigid-body, and is termed the body- fixed frame. The other is fixed in 
inertial space, and is called the inertial reference frame. For motions in which there is 
no translation, a point on the body always remain fixed in inertial space. This fixed 
point can be taken as the common origin of both the frames. 

With this, the attitude of the body may be defined as the relative orientation 
of the two reference frames. We define the The reference attitude of the body as that 
attitude for which these two frames coincide. Other kinematic quantities like angular 
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velocity, acceleration etc. are similarly the relative velocities etc of the bodyfixed 
frame relative to the inertial frame. 

1.1.1 representation of attitude 

The attitude of the body can be represented mathematically by any set of parameters 
describing the relative orientation between these two frames. Based on this, there 
are a few conventions of attitude representation. Most important among these are 
mentioned as follows. 

1. Euler angles : This is one of the most used representations of attitude. It 
represents the attitude in terms of three successive rotations of the body-fixed 
frame with respect to the inertial frame, where the rotations are about the 
body-fixed axes. The angles of these rotations, denoted by 0, 6 and i/), are 
called the ‘Euler angles’. 

2. Quaternions :This method of attitude parametrization is based on ‘Euler’s 
theorem for rigid-bodies’. It represents the attitude in terms of a quadruple of 
real numbers, or to be more precise, a scalar and a 3-Dim vector. Due to its 
suitability in numerical processing and some other desirable properties, it has 
come to be widely used in recent years, especially in spacecraft applications. 

3. Cayley- Rodrigues parameters : This is another 4-parameter represention 
of 3-Dim attitude. Closely related to quaternions and Cayley-Klein parameters 
which also share the fast-numerical analysis advantage like the quaternions. 
They are applied often in other branches like physics wher rotaional symmetry 
is present. 
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1.1.2 Mechanisms for attitude control 

In space applications, various types of mechanisms are used to produce the control 
torques required for executing maneuvers. Two major categories of these are 

• Reaction wheels 

• Retro rockets 

In reaction-wheel mechanisms, the spacecraft has mechanical wheels along the three 
axes. These are called reaction-wheels or momentum-transfer-wheels. By transfering 
the spacecraft’s momentum to these wheels, thereby redistributing the total momen- 
tum along the three axes, a desired attitude of the body can be achieved [3]. 

In the retro-rocket method, small rockets or jets on the spacecraft are fired 
to deliver directional thrusts to reorient the body, these jets are mostly either pulse- 
width-modulated (PWM), pulse- frequency-modulated (PFM), or a combination of 
both (PWPF) [6]. 

1.2 Previous work 

The control of attitude and motion of a spacecraft has been an area of considerable 
investigation for quite a long time. A fair amount of published literature can be 
found, which describes the various approaches regarding attitude representation, and 
mechanisms used. 

In the area of attitude representation [4] gives two numerically efficient al- 
gorithms for determination of spacecraft attitude from vector observations. A basic 
introduction to the use of quaternions for spacecraft attitude representation is given 
in [2j. An algorithm for global attitude estimation in quaternions can be found in [5]. 

The specific problem of stopping a rotating spacecraft while optimizing an 
energy function, has been described in [13]. In [6], it was shown that for attitude- 
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feedback in terms of quaternions, a large angle stable rest to rest maneuver can be 
achieved. A fairly detailed study of the attitude control problem and various control 
laws using quaternions is done in [11]. In [10] a linearization of the attitude error 
dynamics is shown, while [9] gives a control law in which the angular rates need not be 
actually measured, thus requiring no velocity sensors. A globally' assymptodically sta- 
ble control law that takes the spacecraft from any arbitrary orientation to any other 
orientation using quaternion feedback, was suggested in [8]. But there, the gains 
selected were arbitrary, with no optimization procedure suggested. In [7], however, 
a control law achieving a large-angle maneuver through the shortest-angular-path 
Eigen-axis rotation was suggested. This is optimal in path-length, but no method 
or possibility of optimization among such laws with respect to energy or time was 
suggested. 


1.3 Preview of problem considered in this thesis 

In this thesis, we derive control laws for taking a spacecraft from any arbitrary initial 
attitude to any other desired attitude, through an Eigen-axis rotation. The form of 
control law derived is a generalization of the control law suggested in [7]. The form 
suggests a family of possible control laws for achieving the above objective. Possible 
scopes of optimization over this family are also suggested. 

1.4. Organization of this thesis 

In Chapter 2, a brief introduction to quaternions and their use in attitude represen- 
tation and rotational kinematics is given. The Kinematic and Dynamic equations 
governing rotating rigid-bodies are presented in terms of quaternion parameters. The 
concept of Eigen-axis rotation is introduced and its mathematical representation is 
derived. 
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In Chapter 3, a detailed study into the structure of Control laws required for Eigen- 
axis rotation is done. This leads to suggestion of a number of control laws for the 
above purpose. Stability analysis of those Control laws is done, followed by a brief 
robustness analysis. 

Chapter 4 presents the results of Computer simulation of the Close-loop Control sys- 
tems with the developed Control laws along with brief discussions on them. 

In Chapter 5 suggestions regarding possible future work are given. 



Chapter 2 


Quaternion representaion of 
system equations 


In this chapter we present the kinematic and dynamic equations of a rigid-body mov- 
ing with one point fixed, in terms of quaternions. Traditionally, rigid-body kinematics 
is handled through the Euler-angle formulation of attitude. Hence, the kinematic and 
dynamic equations have been given in terms of the Euler angles <p, 6 and "0 along with 
(jj i, the angular velocity components along any of the two frames, body-fixed or iner- 
tial [1]. However these equations involve trigonometric functions of 0, 9 and 0 which 
hinder fast and efficient numerical processing. In addition, they also involve singular- 
ities in attitude representation. Nevertheless, despite these drawbacks, Euler-angles 
have remained one of the most popular formulations. 

The kinematic equations in terms of quaternions as attitude parameters was first 
given in 1958. For a short history of quaternionic attitude representation and kine- 
matic equations, the reader is referred to [7] and [2]. In recent years, quaternions 
have been commonly used as attitude parameters mainly in space applications due 
to some of their important and useful properties. These can be mentioned as follows. 
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1. They have no inherent goemetrial singularities, as do Euler-angles.[ 6] [8] [9] 

[ 11 ] 

2. They have no singularities in their kinematic differential equations, as do the 
Cayley-Rodrigues parameters. [6] 

3. Successive rotations can be represented as successive multiplications by 4 x 4 
quaternion matrices. [2] [6] [9] 

4. They are algebraic quantities, with only products appearing in all their equa- 
tions, without any trigonometric quantities. This makes them suitable for fast 
numerical processing on board spacecrafts. [2] [6] [8] [11] 

The above factors make them ideal for use in space applications. 

2.1 Euler’s theorem and Euler-rotations 

The quaternion approach to representation of attitude and rotations is based on the 
well Euler’s theorem for rigid body motion. The theorem can be stated as follows 

Theorem 1 (Euler’s theorem for rigid-bodies) ' The most general displacement 
of a rigid body , with one point fixed , is a rotation about some axis . ’ 

Let a rigid body execute a single, or a sequence of maneuvers, such that one point on 
the body always remain fixed in an inertial frame. Then the above theorem implies 
that the net resultant transformation from the initial orientation to the final, is always 
equivalent to one single rotation about an axis passing through the fixed point. The 
theorem also implies that for a rigid body moving with one point fixed, the resultant 
attitude after any number of transformations starting from the reference attitude, can 
always be expressed as the result of a single rotational maneuver from the reference to 
the final present attitude. The axis of this rotation is called the Eigen-axis or Euler- 
axis for that attitude. The point that remains fixed obviously lies on this axis. A 
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rotation about the Eigcu-axis, is called an Eigen-axis rotation or Euler rotation. The 
body then be brought back to the initial or reference orientation by the same path in 
reverse. That is, rotating the body about the Eigen-axis would enable it to return to 
the original orientation through a single rotational maneuver, irrespective of the type 
and number of maneuvers that took the body to the final orientation. The physical 
facts mentioned above can also be understood from a mathematical view-point. For 
this, and a detailed discussion of the theorem, the reader is reffered to the Appendix. 
It can be shown that for a given initial or reference attitude and a final attitude, the 
Eigen-axis is unique, irrespective of the intermediate maneuvers. (Appendix) Hence, 
for a given reference attitude, any other attitude can be uniquely specified in terms 
of its Eigen-axis and the angle of the Eigen-axis rotation. Thus these two can serve 
as parametrization of attitude for a rigid-body moving with one point fixed. 

The case is equivalent to the linear case of a point-object moving in 3-Dim space 
starting from a location and reaching a final location after a scries of moves. In 
Fig 2.1 Let Pq be the initial position, and a reference frame {X — Y — Z) fixed with 
Po as the origin. Now the displacement of the body from initial Pq to final Pn can 
be aways be considered equivalent to three consecutive independant displacements 
parallel to the three axes, as shown. In particular they can be reversed to return the 
body to Po irrespective of how it has actually come to the point P„. Hence the lengths 
of these three moves are properties of the location, and can be used to spaecify it — 
being called the co-ordinates of the location. 

2.2 Euler’s theorem and quaternions 

As seen in the previous section, the attitude of a rigid body can be specified in terms 
of an axis, say (ni n 2 n 3 ), and a corresponding rotation about it through an angle, 
say '9', that would take the body from the reference or origin attitutude to the present 
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Figure 2.1: The co-ordinates of a point 


one. 

Based on this, a quadruplet of real numbers given by the following can be used as an 
attitude or orientation measure. 

9 = (9o, 9 ) = ( 90 , 91 . 92 , 93 ) 

where 


and 


go = cos 0/2 

(2.1) 

g = n • sin 0/2 

(2.2) 


n = (ni,n 2 ,n 3 ) 

is a unit vector along the axis of rotation (Eigen-axis). Here g, which is clearly 
collinear with h is called the vector part of the quaternion , and go the scalar part. 
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Equations 2.1 and 2.2 can be expanded into the form 

go = cos 0/2 
Qi = Ui sin 0/2 
ga = ua sin 0/2 
gs = ns sin 0/2 

It can be checked that 

Qo + 9? + ga + 93 = 1 

Thus in the representation of attitude, the quaternions arc of magnitude one, that 
is, they are unit quaternions . It is easy to see that the original or reference ori- 
entation of the body would correspond to the case 0 = 0. with fi being any vector 
whatsoever. Physically, it means the body has not undergone any rotation about any 
axis whatsoever. For this attitude, irrespective of n, the quaternion representation is 
always 

ge = (go, g)c = (i, 0,0,0) 

Various ways of attitude detremination in terms of quaternions, using an inertial ref- 
erence framc,can be found in [4], [2], and [12]. In short, the relation between the 
quaternion representation and the Direction-cosine and Euler-angle representations, 
can be understood through the following observations. The orthogonal transforma- 
tion matrix or the rotation matrix, corresponding to a rotation about the axis n 
through an angle 0 is given by 

g? + go-2 9192-9390 9 i 93 + g2go 

9192 + 9390 92 + 90 “ 2 9293 — 9i9o 

9391 - 9290 9392 + 9i9o 91 + 9o “ 5 


A = 2 


(2.3) 
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where go, gi, ga, <73 are related to 9 and n by the Eqs 2.1 and 2.2. The above can be 
derived from a basic geometric consideration of the rotation (Appendix) and is related 
to the well-known ‘rotation-formula’ and the ‘Rodrigues formula ’ in mechanics. 
Conversely, if the rotation matrix A is given, say as the Direction-cosine matrix itself 
or in terms of Euler-angles , then the vector h is easily recognised as the Eigen- vector 
corresponding to the Eigen-value 1 of the matrix and the angle 9 is found by noting 
that the other two Eigen- values are e*" and ( Appendix ). Thus, by considering 
Eqs 2.1 and 2.2, g may be found . 

A much more direct way of calculating g, form the rotation matrix A is sugested by 
Eq 2.3. Looking at Eq 2.3, we can observe that 


1 + trace[A] = 4 ql 


By using the constraint 


(Jo + 7? + g 2 + g| = 1 


we have 


go = -y^l -f trace[A] 


Further, it can be seen that 


A — AF — 4 go Q > 


where Qx is the cross-product matrix of g , given by 

0 -73 72 

= qx = 53 0 -gi 

-72 7i 0 

Thus, from the above, and using Eq 2.4, we have 

1 1 


Qx — 


2 a/I -I- tracelA] 


{A - A^) 


(2.4) 


(2.5) 
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Hence, both qo and q are found. Generally, both q and —q may represent A as 
mentioned in [9], but the kinematic equations considered here (Eq 2.7), as in [9], are 
consistent with q representing A. A numerically efficient algorithm for computing q 
from the rotation matrix can be found in [12]. 

The above discussion shows that the numerical and control processings can be done 
in terms of quaternions, even if the attitude is measured actually in terms of other 
parameters like Euler-angles. In these cases, the above relation*, or the algorithms in 
[9] may be used to transform the attitude representations to quaternions. 

2.3 The Control Problem 

In most control applcations, the control laws are designed to take the body to an 
orientation for which the quaternion representation is (1,0, 0,0). As seen in the 
previous section, this is a representation of the reference attitude of the body. When 
this is the target value of the quaternion q = (qo,q), the value of q at an instant 
is itself an error measure. Hence, the quaternion qc = (1,0, 0,0, ) is taken as the 
‘origin’ or ‘zero’ attitude, and the control problem formulated as a regulator problem. 
Control laws are developed which take the spacecraft from a non-zero orientation , 
represented by 

9(0) = (ffl(O) , 9(0)) ; 9(0) 5^ 0 

to the origin or ’zero’ orientation 

qe = ( 1 , 0 , 0 , 0 ) 

However this does not restrict the generality of the attainable target attitude. In 
other words, this does not imply that the control laws will take the spacecraft always 
to a fixed particular orientation in 3-Dimentional space. This is because any targetted 
orientation of the body can be taken as the reference attitude corresponding to the 
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origin qe> a-nd the present orientation g(0) is expressed with reference to this, thus 
having non-zero f7(0). Hence the problem taking the body from any given present 
orientation to any other target orientation, is translated to one in which it is taken 
from an initial attitude ^(0) with non-zero g(0) to (1,0, 0,0). 

This formulation is an often-used one and has been used in [6], [7] and [8]. A detailed 
discussion regarding possible other attitude error measures using quaternions can be 
found in [11] . ^ 

2.4 Kinematic and Dynamical Equations 

In terms of conventional Euler angles ^ , d and the kinematic equations can be 
written as follows. 

(p = (sin i/j / sin 0) u)i + (cos V; / sin (?) U 2 

6 = (cosi/j) u)i — (sini/’) W 2 

■ip =1 — (sin V-’ cot 0) tui — (cos ip cot 9) + u)z 

Here uji , and uz are the angular velocities about the three body-axes. It 
should be noted that similar equations can be written where cuj are components of 
the angular velocity in inertial frame. Also, various conventions of taking the Euler- 
angles give rise to various similar but different equations [Ij. The dynamic equations 
relating the angular-velocities and the external torques are given by the well-known 
Euler equations of rigid-body motion. For a general case when the body-fixed axes 
may not coincide with the principal axes, they are written as 

Ju) = —io X Jw •+• r (2.6) 

where Jsxs is the inertia matrix, u the angular velocity vector and r the torque 
vector. A detailed discussion and derivation of the above equations can be found in 
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any standard treatise on rigid-body mechanics . 

Next, we consider the kinematic equations of rigid-body motion with one point 
fixed, where the attitude parameters are quaternions. The are stated as follows 

q = \F[q).Q (2.7) 

where 

<7 = [ (7o <7i 72 qz V 
Q = [ OJi U2 Us Y 

and F{q) is the matrix given by 

-ql -q 2 -qs 

Qo 93 -92 

-Qz 9o 9i 
92 -9i 9o 

It may be observed that the quaternion-kinematic equation 2.7 involves only alge- 
braic terms, while those for Euler-anglcs involve trigonometric terms, illustrating the 
numerical advantage of quaternions. 

Further, the matrix F{q) can be seen to satisfy the following important relations 

■ Hi) = 4x3 ; F^{q) • f = 0 

Using the above, and Eq 2.7, we can obtain 

9^ 9 = ■ F{q) • (D = 0 • 0) = 0 

which implies 


m 


-T 

■9 


Qohxz — Qx 


q^ q = Constant 




2. Quaternion rcpresenUiion of system equations 


17 


This is consistent with the constraint on the norm of q, always being 1. In other 
words, the solutions of Eq 2.7 are all unit-quaternions, representing attitude. 

Eq 2.7 can be written separately for the scalar and vector parts, as 


1 -T- 
--9 O) 

(2.8) 

-qoQ -Q^u 

(2.9) 

2 

-qoQ-qxu! 

(2.10) 


The dynamical equations for rigid rotating bodies are the same as in Euler-angle 
formulations, that is same as Eq 2.6. Hence, the Eqs 2.8, 2.10 and 2.6 are the 
equations describing the spacecraft motion. Then, considering the torque vector r as 
the input vector the control system equations in terms of quaternions are 

Qo = -\f^ (2-11) 

q = -f (2.12) 

J i) = —ClJu) -i- u (2.13) 


where 0 is the cross-product matrix of cn, given by 


n = 


0 - 0)3 UJ2 

U )3 0 — 0)1 

- 0)2 0)1 0 


The equations no 2.11, 2. 12, and 2.13 give our system state-space equations , with 


. ?o Ql <73 ^2 


as the states and 


’ Ui U2 Uz 


as the controls. 
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2.5 Euler-rotation and its mathematical 
expression 

We expalained above how a manuever of a rigid-body with one point fixed from 
one attitude to any other can be achieved by a single rotation of the body about 
an axis. The axis is called the Eigen-axis or Euler axis, and given such an axis h, a 
rotation of the body about this axis is called an Eigen-axis rotation or Euler-rotation. 
Morever, we can take the target attitude as the reference attitude with the quaternion 
representation (1,0, 0,0). The above formulation is used for developing control laws 
for changing the attitude of the attitude while the initial attitude is represented by 
an arbitrary quaternion 

g(o) = ( <70 ( 0 ) . m 

O’ I 


o 



Figure 2.2: Eigen-axis rotation 

Referring to the Fig 2.2, let the initial attitude of the body be depicted by a 
certain position vector OP' of a point P' while the target attitude has the same 
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position vector rotated to OP through an angle ^o- 00' is the Eigen-axis, with unit 
vector n and the rotation trajectory is as shown. Recall that the vector ^(0) is related 
to n and the angle d of the Eigen-axis rotation through the Eq 2.1 and 2.2. 

As the body traces the shown Eigen-axis trajectory to the target orientation OP , 
it is clear from the figure that all intermediate orientations also have the same Eigen- 
axis n, while they have coiitiniously decreasing 6 from 9o to 0. Hence , The attitude 
of the body undergoing Eigen-axis rotation will be such that the instanteneous vector 
part of the quaternion is given by 

q{t) = n-sin9{t)/2 

which always remains along h . The direction of the angular velocity vector, will also 
evidently remain along the direction of n or g throughout, the magnitude starting 
from zero and returning to zero at the target point . 

Hence , we can write , for the body to execute Eigen-axis rotation, the solution to 
the Kinematic and Dynamic equations 2.11,2.12 and 2.13should have the form 

q = Cg{t) ■ q{0) 
uJ = Cu (t) • g(0) 

with 

^.(O) = 1 ; = 0 

a(0) = 0 ; a(t/) = 0 

where tf is the final time. 

The objective of the control law for an Eigen-axis rotation is to achieve the above 
solutions of the control system equations 2.11,2.12 and 2.13 once the control law has 
been substituted into the eq 2.13. A detailed analysis of such control laws is carried 
out in the next chapter. 




Chapter 3 


Closed-loop Control laws for 
Eigen- axis Rotation 


In this chapter we study the structure of a special form of control law that takes a 
spacecraft from any initial orientation to any given target orientation, through an 
Eigen-axis rotation as described in the previous chapter. A broad outline of the 
formulation of the problem was also given in that chapter. Based on this, we develop 
few control laws that are found to be generalizations of a control law given in [7] 
for the same objective. Global asymptotic stability for these control laws are shown, 
and where required, the conditions on these control laws needed to guarantee global 
asymptotic: .stability arc also cleterininccl. 

3.1 Developing the Control law 

In [7], it was shown that a large angle rest-to-rest manuever about the Eigen-axis can 
be achieved by a quaternion feedback control law, with proper selection of feedback 
gain matrices of the quaternion feedback regulator. The control law suggested fir this 
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was of the form 


u = n./w — Du) — Kq 


where the gain matrices K and D were given by 

K = kJ ; D = dJ 

k and d being two scalar constants. In that work, however no analysis leading to the 
suggested form was given. Further, it appeared that the restriction of the possible 
gain matrices to constant multiples of the inertia matrix was rather severe. In this 
chapter we try to cirrivc at a more general form of the control law required to achieve 
an Eigen-axis rotation, starting the analysis from the required form of the Eigen- 
axis trajectory. The more general form that we try to develop consists of general 
time-dependent gains. 

3.1.1 The Control Problem 

The quaternion kinematical equations are 

qo = -(l/2)r7^’a; 

^ = (1/2) nf -I- {1/2) QoQ 

where ql + = 1 ; 

and n represents the cross-product matrix of the angular velocity vector uj 

0 — W 3 U)2 

^ = UI 3 0 —u>i 

—UJ2 Wi 0 

The dynamical Euler’s equation arc given by 



Ju! = —QJu! + u 
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where J is the 3x3 Inertia matrix and u the external torques. With this the control- 
system equations are 


go = -(1/2) 


(3.1) 


^ = (1/2) Qg -h (1/2) go w 


(3.2) 


Ju = —QJu -f u 


(3.3) 


with 


{Q0,q,djV = (?0,<?l)?2,g3,‘^ l,w 2,W 3 )^ 


as the state vector, and 


U = (til,U2,'U3)^ 


as the input vector. 

\ 

The problem considered is to find input u so that the solution to the above system 
is of the form 


■ ?(^) 

(3.4) 

(7a, {t) • g(0) 

(3.5) 


where g(0) is the initial value of the vector q and Cg(t) and Cb (t) are scalar functions 
of time, with 


(7,(0) = 1 ; C,{tf) = 0 

a (0) = 0 : C^ {tf) = 0 


where t/ is the final time. 
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3.1.2 A special form of the control Icfw 

We choose a general feedback law involving the quaternion vector and the angular 
velocity, given by 

u = fiJw — Mill — Nq (3.6) 

Where M and N are 3x3 time-varying matrices . We need to find suitable M and N 
so that with the above input the solutions to the close-loop system equations describe 
an Eigen-axi.s rotation. Substituting Eq 3.C into Eq 3.3, it gives 

J-^Nq 

or 

cD = —AG) — Bq (3.7) 

where 

A = J-^M ; B = J~^N 

Hence, Eq 3.1, 3.2 and 3.7 are the state equations describing the close-loop control 
system which should have solutions of the form 

q = Cq{t) ■ q{Q) 

dJ = a(t)-g(0) 

Hence, substituting these in Eqs 3.2 and 3.7, the necessary conditions for Eigen-Axis 
Rotation become the existence and Cq such that 

Cq q{0) = ICqC^ q{0) x g(0) + yoC^q{0) 

■ Gq = 

= - O’C? 

with 


Q = II m f 
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and 

a = -AC^-BCg 

Thus in order to have an Eigen-axis rotation , the following equations should be 
satisfied 

c, = 590 a = ic„^l - Q>C| (3.8) 

CJ = -AC^ - BC, (3.9) 

where 1 is the 3x3 identity matrix. We need to ascertain the conditions on A and 
B so that a scalar solution to Eq 3.9 exists . 

3.1.3 Structure of matrices A and B 

To investigate the structure of matrices A and B, we breakup Eq 3.9 into two parts, 
one involving the diagonal elements and the other involving the off-diagonal elements. 
We write them separately as 

C'u, = -an a - ki Cg for i = 1, 2, 3 (3.10) 

and 

aij + kj Cg = 0 for i = 1, 2,3 ^ i (3-11) 

We analyse the two separately fis below. 

The diagonal elements : 

The three equations 

Cfjj — ttji Cu bn Cq 

will have a unique solution if and only if 

= UuCk( + biiCq — OllOu; + = A (3.12) 
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where A is an arbitrary variable difrereiitiablc with respect to time. Thus the above, 
is a condition which must be obeyed by the diagonal elements of matrices A and B 
for achieving Eigen- axis rotation. Now, for developing such control laws, we observe 
that Eq 3.12 can be satisfied by the diagonal elements in any of the two following 
cases. 

(a) 

flu = 0.22 — <^33 = 

/)] 1 = 1)22 = ^33 = 

where Od and bd are two arbitrary scalar variables, differentiable with respect to time. 
For this, Eq 3.10 becomes 

Cu ~ — Od C(j} — bd Cq (3.10 b) 

(b) 

bii — {A 0>ii\/Cq 

where A is any time-differentiable function. For this case, Eq3.10 becomes 

=-A (3.10 c) 

The off-diagonal elements : 

The six equations in Eq 3.11 

\ 

o-ij Cui + C, = 0 for i ^ y ; i, j = 1, 2, 3 

can be satisfied only through the following two cases. 

(a) Oij — bij — 0 for i ^ j , — 1, 2, 3 

(b) bij =■-%-■ a.ij i,i = 1, 2, 3 

Any of the above will ensure satisfaction of Eq 3.11, which is required for an Eigen- 
axis rotation. But it should be noted that Eq 3.11 implies that the contribution of the 



3. Closed-loop Control laws for Eigen-axis Rotation 


26 


ofF-cliagonal elements to the right-hand expression of Eq 3.9 is zero. Thus in either 
of two cases above, the off diagonal elements Cy and bij for i ^ j do not in any 
way contribute to the dynamics of the system. Therefore it is reasonable to consider, 
for developing control laws, only the case (a) for which the off-diagonal elements are 
simply taken as zero. 

Hence, finally combining the results of the analysis of the diagonal and off-diagonal 
parts as above, we get two sufficient sets of conditions on the elements of A and B 
for achieving Eigen-axis rotation. They are given below. 

( 1 ) 


dn = ^22 

= azz = 

dd 

= ^22 

= 633 = 

bd 

(Xij 

= % = 

0 

bij = {A — 

dijCfjj }/ C /(7 


dij = bij 


= 0 


For this thesis, we will consider only the simpler form of the control law arising 
from the conditions of set (1) above. The matrices A and B that follow from these 
conditions are 

A = aj = J-^M 
B =■ bj = J-^N 

The above gives a condition on M and N , so that the body undergoes an Eigen- 
axis rotation . Changing a notation , writing two generally time- varying scalar func- 
tions ip[-) and (i>{-) in place of ad and bd , we finally arrive at 


M = i}{-) J 
N = <!>{■) J 
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Thus the result is proved and we can write the control laws that achieve Eigen-axis 
rotation as 


u = Q. JQ — i/jJoi — 4>Jq (3.13) 


A block diagram of the resulting Control System is shown in Fig 3.1 



Figure 3.1: The Control System 


3.1.4 Remarks 

It may be mentioned that the first term was introduced to cancel the cross-coupling 
torque. In most cases this is found to be generally much smaller than the other 
torques involved [7]. As a result, even the exclusion of this term does not affect system 
response significantly. This fact was verified in our simulations of the above control 
laws. The reason behind the above is that in most cases, the principal moments of 
inertia of the spacecraft are considerably larger than the cross moments of inertia, 
as a result of their being generally symetrical . In addition, the body fixed axes are 
taken to coincide with the body’s primary axes. In either case the inertia matrix J 
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is nearly or completely diagonal, making the vectors Q and JcD nearly or completely 
parallel. As a consequence the term Q Ju, which is w x JQ, is small or zero. 
Therefore, this term is required only in those cases where this cross-coupling torque 
is significantly large. 

We thus conclude that the control law of Eq 3.6 need not be restricted to those with 
gain matrices being a constant multiple of the inertia matrix. We can have a family of 
control laws with different time-varying cj) and i/*, They can be purely time-functions 
or may be also be taken as scaler functions of the state {q,0 ). The control law 
suggested in [7] is evidently a special case fo this, with both (j) and -ij} as constants. 
This leaves us with a scope of optimization over a number of control laws, with respect 
to some objective function. This also requires a study of the conditions on (j) and if), 
for the system to be stable. 

3.2 Stability analysis 

In this section, we determine conditions on time-varying functions 4>{t) and i/)(t) so 
that the close-loop system has the i)oint (</o = 1,9 = 0,d7 = 0) as a stable equi- 
librium. For this we make use of Lyapunov’s second or direct method for a general 
time-varying system. The methosd is based on the Lyapunov theorem for such sys- 
tems, a non-rigorous version of which can be stated as follows. 

Theorem 2 ( Lyapunov theorem) A general time-varuing system described by the 
equation 

^-f{x,t) 

has the point as a sto.ble equilibrium point, if in a neighbourhood of x^ there 
exists a continious scalar function V [x, t) such that the folowing conditions are sat- 
isfied. 
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1. 0 = 0 for alU 

2. y(a;, t) > 0 for all x 7^ Xa in the neighbourhood Nx^ for all t 

3 . V{x, t) <0 along the trajectories of the system in Nx^ for all t 

If the above is true for the entire state-space, i.e. if Nx^ is the entire state-space, then 
Xf. is globally stable. If condition ( 3 ) is satisfied as a strict inequality, i.e. V(x,t) < 0 
for all points other than the equilibrium, then the point Xg is asymptotically stable. 

In our problem the system equations are given by Eq 3. 1,3.2 and 3 . 3 . After 
substitution of the control law given by Eq 3 . 13 , we get the close-loop system as 
below 


9o = -(l/2)g^cj 

(3.14) 

g" =(l/2)ag--f (1/2) go CD 

(3.15) 

to = — — (f)q 

(3.16) 


The equilibrium point is obviously the target attitude 

9o = 1 ; ?! = 92 = 93 = 0 ; = 0 ; 

We can see that a simple positive scalar function of the instanteneous state that 
vanishes only at the equilibrium point can be given by 

e = (l/ 2 )d;^w + (1 - 9 o)^ + 9 i + 92 + 93 
= {l/2)u)^u) -h 2(1 — go) 

Suggested by this, we consider the following two similar Liapunov functions of the 
form 

Vi = I ui -f 2(^(1 — go) 

^ “ 5 i + 2(1 — go) 


and 
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with the restriction 


(^(t) > 0. 

Hence, we have 

) ^2 — 0 at the Eq. Point, for all t 

^i,y 2 > 0 for all other points, for all i 

Differentiating Vi, and using the system equations 3.15 and 3.16, we have 

Vi = oj'^Q — 24>qQ — 2(^(1 — (Jo) 

= gF{-'>Puj - ^q] - 2(j){-{l /2)u}^ q] + 2^(1 - go) 

= — -0 {uFq) — <f>Q'^q + (f)U3^q + 20(1 — gO) 

= -0 {ui^u) + 20(1 - go) 

Looking at the above expression , and noting that go is always > 0 . We conclude 
that the equilibrium point is globally asymptotically stable if the following set of 
conditions is satisfied 

0 > 0 

0 > 0 

and 0 < 0 

Thus , the above is a set of sufficient conditions for global asymptotic stability, 
emerging from the first Liapunov Function Vi . Now, differentiating V 2 , we have 

= 2 ' ~ ~ 2go 

= - r -2{-|dJ^g} 

= ^{-0w - 0g} - ■ cJ + cJ ^g 

= ^ • cJ — (J Tq — ^ ■ -^uJ '^dJ +6J^g 
~ — 5 (200 + 0)aJ '^uJ 

From this we note that the equilibrium point is globally asymptotically stable if 

(200 + 0 ) > 0 
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A sufficient condition that ensures this, is 

(j) , 'Ip >0 with p > 0 

Thus another set of sufficient conditions, emerging from the second Liapunov Func- 
tion is 

<^ > 0 
p > 0 
and P > 0 

Combining tlicsc two sots of conditions, oach of which is sufficient to guarantee 
global stability, wc sec that as long as p > 0 and p > 0, both the cases of > 0 
and p < 0 ensures global asymptotic stability. Hence, the only required conditions 
sufficient to guarantee stability, are 

p > 0 
p > 0 

The control law proposed in [7] is clearly a special case of the above, with p and p 
as positive scalars. 

3.3 Other control laws 

In the previous sections we looked into the form of the feedback controls that ensure 
that the spacecraft executes an Eigen-axis rotation. The form of control laws we 
arrived at were 

u = flJcj — pJu) — pjq 

where p{-) and p{-) were scalar variables. We also obtained sufficient conditions on 
controls p and p that ensure global asymptotic stability, in the case when they are 
time- varying. In general, however, p and p can be scalar functions of the state- 
variables. In this section, a few more control laws of the above form are suggested 
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and tested. They are special Cciscs of the above general form of control law, with (f), 

Ip, or both S])ociiicd. Global asymptotic stability of the desired target point for these 
control laws is also demonstrated. 

It may be mentioned that the purpose of suggesting the forms of ^ and ip in these 
control laws is primarily to illustrate the general nature of the above functions while 
achieving the Eigen-axis rotation. It illustrates the possibility of taking a wide choice 
of functions of the state variables as (p and ip for achieving Eigen-axis rotation to 
the target attitude. The suggested control laws are examples from a possibly large 
number of control-laws that lead to an Eigen-axis rotation from any initial orientation 
to any target orientation. 

In general, the system equations for the rigid-spacecraft rotations are of 7-dimentions 
with actually 6 independant equations, the 1st of the quaternion kinematic equations 
being derivable from the other three, using 

fQ + 9o = 1 

But with control laws of the above form the body undergoes an Eigen axis trajectory, 
along which q and Q are scalar multiples of the initial g(0) as given in Eq 3.4 and 
Eq 3.5. By substituting those in the close-loop Equations 3.14 3.15 and 3.16, we get 
the following equations 


= lc„yi - Q 2 C| 

(3.17) 

= -0 c„ - c, 

(3.18) 


In this case the system equations get reduced to those of a two-dimentional system, 
with Cq and as the state variables. With (p and V-' remaining scalar, the above 
two-dimentional system is always equivalent to the 7-dimentional original system. 
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with the scalars Cq and being proportional to the magnitudes of the q and dJ , 
whose directions have been determined by the initial attitude q{0). We can consider 
the above equations to be a two dimentiomnal system, with Cq and Ccj as the state 
variables, and 0(-) and ^/'(•) as the controls. Within this reduced system, and. 
ip{-) can be pure functions of time, in which case the reduced system is an ‘open-loop’ 
one, while the original system can be said to be a combined ‘open-loop - closed-loop’ 
one. On the other hand, if 0(-) and ■?/)(•) are functions of q and/or d), then the reduced 
system becomes a closed-loop controlled one, while the original system can be looked 
upon as one with a highly non-linear feedback. 

With these ob.servations, we proceed to suggest and analyse a few control laws. 
3.3.0. 1 Control law 2 

A control law is suggested by taking (f> and tp as below. 

0(0 = </ o /2 : ' 0(0 > 0 

With the abov('. control law the general close-loop system becomes, 


(jo = -{1/2) if u 

(3.19) 

q = {1/2) Clq -t- {l/2)qoui 

(3.20) 

Ji) = -0 Ju) — ^ (Jo ? 

(3.21) 

. / - 1 - 
w = -Ip u - -qoq 

(3.22) 


The equilibrium point is obviously, 


q = u} = 0 
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The stability ol the equilibrium point is studied with the following Lyapunov function. 

V = + ifq 

= 0 atEq.(g = 0 cJ =0 
> 0 for all other points 

Now differentiating this, and using Eqs 3.20 and 3.21, we have 
V = tJ '^(j + q'^q 

= i/'dJ -5909} + } 

= - i)Ld - \q^uj 4 - 0 + \qnfu 
= — '0 dJ dJ 

< 0 for all points other than the Equilibrium 

Thus the system is globally assymptodically stable, with respect to the point 
5 = (D = 0 Simulation results show that with suitable choice of or the 
above control law achieves Eigen-axis rotation to the target point reasonably fast, with 
amount of energy spent being also reasonable - being comparable to that spent for 
low constant values of 0 and The only considerable drawback of this form of 0 is 
that for initial displacement of 0 = 180° exactly , in which case , qo = cos 90 = 0 
, along with d) = 0 , the control torque becomes zero, and the system fails to move. 
In other words, 

<7o = 0 

d) = 0 

is another equbilibrium point, although unstable . Hence, even for angles very close 
to 180° , the body returns to the target equilibrium point satisfactorily. 

3. 3. 0.2 Control law 3 


In this control law, we take 0 and "0 as below. 
0(.) = (1 -f- 7d)'^’(7) ; 0(.) > 0 
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Another lorin of the {mictions that was tried was the above, where 7 is a +ue scalar. 
We consider the lyapunov function 

+ 2(1 — 50) 

= 0 atEq. {5 = 0 uJ =0) 

> 0 for all other points 

The equality occuring only at cJ = g = 0 . 

Differentiating , and using the ( ) , we have 

V = u ~ 2 (fo 

= - {I + ju)^q)q} — 2{-\dj q} 

= —'ipcJ ^ dJ — Q'^q — fj)^ + tJ q 
= —xpcj’^uJ — 7 (d 7 ^g)^ 

< 0 for all points other than the Equilibrium 

Thus the system is stable, globally and assymptodically, with respect to the point 

Q = 0 \ q = 0 ] go = I 

at the same time achieving Eigen-axis rotation, which is verified by simulations. 

3.4 Robustness 

The robustness of the above suggested control law, was checked by simulations, 
tlirough variations in the inertia matrix used for feedback. Sufficient robustness was 
observed in all of them, for considerable variations. The nature of robustness for one 
of the control laws can be seen by the following analysis. Considering the control law 
with 4 > = (/0/2, wc! have 


u = ClJ'u — — 4>J'q 
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where J' is the inertia matrix with some error. Considering the same Lyapunov 
function ixs in control law 2, 

T / ^ — T ^ — T 

y = -w dJ + 

we have 

V = (J ‘\j q'^'q 

= w • JuJ + q^'q 
= dj +QJ'cJ -ipj'uj -\q^J'q} 

+f/'^ + \qQ(jj } 

= uJ J~^Cl {J' — J) (J — 'ipuJ {J~^ J')dJ — ^qodo J~^ J' q 

+ 0 + ^Qoq'^dj 

= - J)u - - igodJ ( J- V' - J) g 

If the error in measuring J, i.e. the difference between J and J', is not unusually 
large, then both (J' - ./) and (J~‘ J' — I) arc sullicicntly small, and J~^J’ is near 
about / 3 x 3 - At least, (J' - J) and - 1) are considerably smaller than J“V'. 

Hence the middle term in the above expression dominates. Morever, for not large 
errors between J and J', (J“^ J' - I) is +ve definite, making the expression for V - 
ue, thereby guaranteeing stability. Even for considerable errors, simulations often show 
stability of the system. However, it was seen that although stability of the equilibrium 
point remained unaffected, considerable errors in the inertia matrix, mainly in the 
diagonal elements, led to deviations from the Eigen-axis trajectory. 

3.5 Conclusion 

The analysis in this chapter led to a particular form of feedback laws that are capable 
of taking a spacecraft from any attitude to any other desired target attitude through 
an Eigen-axis rotation. Within this form, the possibility of having several particular 
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control laws wore also illustrated, through three specific examples. That these control 
laws achieve the above objective have been clear from the mathematical analysis. 
Their voriheation with a computer simulation of the control system and its response 
for the control laws are done in the next chapter. 



Chapter 4 


Results and Discussions 


The main objective of this thesis was to investigate the nature of control laws that 
can achieve an Eigen-axis rotation, from any arbitrary orientation to any other target 
orientation. The emphasys has been given towards highlighting the possibilty of 
generalising such control laws mentioned in previous literature. The control laws 
studied in this thesis were simulated on a MATLAB package suitable for efficient 
numerical procosvsing. The results of these are presented in this chapter. 


4.1 The control laws 

It may be recalled that Eq , , and , gave the control system equations, with u being 
substituted by the proposed control law of the form 

u = flJu) - tp{-)u - 

The control laws analysed may be recalled here. 

1. <l){t) , •0(0 > 0 

2. 0 = go/2 ; 0(-) > 0 

3 . 0 = (1 + 'fuP'q) ; 0 (-) > 0 
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4.2 Response parameters 

Simulations were carried out for the derived control laws, for a wide range of initial 
attitudes given in ternns of their Eigen-axis fi and the angle of Eigen-axis rotation 0. 
In each ol the above control laws, a few different simple forms of •;/>(•) were considered. 
The feature that was common to all the responses was that all of them followed the 
Eigen-axis rof.ation, its expected. This was verified by checking the value of the 
cross-product between < 7 ( 0 ) which is coincident with the Eigen-axis for the initial 
attitude, and the iustanbuiious (inatcrnion vector (]. For Eigen-axis rotation, these 
two vectors have to collinear throughout and hence their cross-product have to zero 
for all instants, which was verified. Apart from this, the other response parameters 
that were observed were the settling time, and nature of response in terms in terms of 
oscillatory or non-oscillatory. For a comparative observation, in some cases a measure 
of the total Eiuirgy involved in the manuever was also considered. This was taken as 
the integral of Ihe body’s Kinetic Energy through the manuever given by 

^ ■ (It 

2./0 

where i.s the body’s Kinetic energy and tj the final time. 

4.3 General observations 

The control law do not impose any condition on ' 0 (-) except that it should be positive. 
In the simulation.s ?/;(•) was taken as function of time alone, which is nevertheless not 
necessary in control laws ( 2 ) and (3). It was observed that for cases where 1 }} rises faster 
compared to <f), the angular rate feedback tends to dominate the feedback after the 
first few seconds. Hence towards the end, with little attitude feedback in terms of the 
quaternion vector but high rate feedback, the body tends to come to a stop fast even 
before reaching the desired attitude (1,0, 0,0). However, if is taken as a positive 
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but decreasing iunction, tJie problem does not seem to arise. Another important fact 
that should be noted is that (j) and V-* leaves sufficient scope for shaping the trajectory. 
For various ty])es of these functions chosen, the nature of the response may vary over 
a wide range. Simulations showed a number of cases in which the response was highly 
o.scillatui'y although allways ixunaining on the Eigen-axis rotation trajectory. 


4.4 Results 

4.4.1 Control law 1 

This control law 

u = njcD — (l){t)q — il}{t)uj 

contains both ^ and as +ve functions of time. The following two typical cases are 
presented. 


(a) ^{t) = 0.5t ; = t 

The above control law was tested with a number of initial attitudes of the body. The 
response in one such case is presented in Figs 4.1 and 4.2. 

The initial attitude in this case was such that it could be reached from the reference 
orientation through a rotation of 0 = 70° about the line n with direction ratios ( 2, 
-4, 3 ). This arbitrarily chosen attitude corresponds to the quaternion ( 0.8192 0.2130 
-0.4260 0.3195 ). The response was smooth, with settling time around 25 sec. The 
cross-product of q and g(0) is shown along with the quaternions in figure 4.1 and was 
found to be zero throughout. 

(b) ^ = t ; 'tp — 0.5t In this case the initial attitude was taken 
as that reached by a rotation of 100° from the reference attitude about the line ( 6 
3 -80) which corresponded to a quaternion of ( 0.6428 0.4402 0.2201 -.5870 ). The 


w1 .w2.w3 ( rad/s ) qO.ql .q2,q3 & qXq(0) 






w1 ,w2,w3 ( rad/s ) qO.q 1 .q2.Q3 A qXq(O) 
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response is shown in Figs 4.3 and 4.4. 



Figure 4.3: Control law 1(b) attitude 



Figure 4.4: Control law 1(b) ang.vel 
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4.4.2 Control law 2 

The control law is 

<^ = ; V' > 0 

This control law luis i) <is a 4-ue function. We present here two simple cases. 

(a) ij) = 0.9 

This is a ca.s<' of constant t/* with ^ as a function of the quaternion. The initial 
condition cosichu’cd corresponded to the quaternion ( 0.9239 -0.3146 0.2002 0.0858 ) 
The response is shown in Fig 4.5 and 4.6. As can be seen, the response is smooth 
non-oscillatory with settling time about 15 sec. The cross-product of q{0) and q is 
again shown along with the quternion responses in Fig 4.5, and is found to be zero 
throughout. 



Figure 4.5: Control law 2(a) attitude 
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Figure 4.6: Control law 2 (a) ang.vel 
(b) t/) = 30/(f + l) 

This is a avsc whcn'o xj> is a clecrcixsiug function of time. The initial attitude cor- 
rcspoiulotl to ( 0.7132 0.4047 0.5665 -0.0809 ). For this, the response is shown in 
Figs 4.7 and 4 . 8 . The response was found to be smooth non-oscillatory. The time 
taken for sc'ttling was around 35 sec. 

4.4.3 Control law 3 

The control law is 

(j) = 1 -ir'yej'^q ; ^(O-O 

where 7 is any positive scalar constant or function. can also be any positive 
scalar function. Evidently, there are a large number of possibilities. We consider the 
following two forms, whose responses to arbitrary initial conditions are presented. 



w1 .w2.w3 ( facias ) q0.q1.<|2.q 3 & qXqCO) 
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Figure 4.7: Control law 2(b) attitude 



Figure 4‘.8; Control law 2(b) ang.vel 
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(a) <l> — I + 0.5 q ; ijj = t 

The initial attiliuie taken for the above is (0.6428 0.6221 0.4397 -0.0733 ). The 
response i.s .shown in Fig.s 4.9 and 4.10. 



Figure 4.9: Control law 3(a) attitude 

Thi.s is a t:}us(! wlicrc towards the later stages, the angular rate feedback dominates 
the control law. Hence, it can be seen that towards the end , the approach of the 
body to ( 1, 0,0,0) slows down. The body almost comes to a stop at attitude slightly 
away from the de.sired (1, 0,0,0,). The response, however, is smooth and the body as 
usual trace.s the Eigen-axis rotation trajectory. 

(b) <l>^l+t-uFq ; V'(-) = l + 50/(t+l) 

This is a case in which tl) has been taken as a decreasing function. The response is 
shown in Fig 4.11 and 4.12. 

The problem of premature stopping is not seen here. The initial attitude arbi- 
trarily chosen corresponded to a case of n having the direction cosines (39 -12 ) 
and e being 175". This implied cfO) =(0.4362 0.1959 0.5878 -0.78.37 V ThP 
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Time t (sec) 


Figure 4.10: Control law 3(a) ang.vel 



Figure 4.11: Control law 3(b) attitude 
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Figure 4.12: Control law 3(b) ang.vel 

was saiuotli and l.iio spacecraft settled comi)Iotely by 52 sec at the desired attitude 
following Eigen-axis rotation. 

4.4.4 A comparative view 

For a rough comparison between the control laws, wc considered the response of them 
for an arbitrarily chosen attitude corresponding to 0 = 120° and h having direction 
cosines ( 17 12 -2 ). For this ^(0) =(0.5000 0.7043 0.4971 -0.0829) Below we give two 
cases, one each for Control laws 2 and 3. It should be noted that the conclusions 
reached may not be general, and may depend on the initial attitude. 

Control law 2(c) : 

<^=§9o ; '(jj = l5/{t + 0.1) 

The responses for this are shown in Figs 4.13 and 4.14. The Kinetic Energy of the 
spacecraft during the maneuver is shown in Fig 4.15. 




w1.w2.w3 < md^s ) Q0,ql.q2.q3 S qXq(O) 
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Figure 4.13: Control law 2(c) attitude 



Figure 4.14: Control law 2(c) ang.vel 
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Figiirc 4.15: The kinetic energy for 2(c) 

The intogrcai of the kinetic energy over the entire maneuver given by | 
was found to bo 700 J. In general, it was found that if the numerator of tp which is 
15 here is incrciisod gradually this total kinetic energy over the maneuver decreases. 
However, the settling time gets increased. 

Control law 3(c) : 

<^ = 1 + ; i> = 10 

The rcsiionsos for this is shown in Fig 4.16 and Fig 4.17. The variation of the Kinetic 
energy of the spacecraft through the maneuver is shown in Fig 4.18. The integral of 
thge kinetic energy for this case was found to be approximately 500 J. 



wl .w2,w3 ( rad/s ) qO q | ^q2.Q3 ^ qXq{0) 
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Figure 4.18: The kinetic energy for 2(c) 

4.5 Conclusion 

In this chapter the responses of the control system for the suggested control laws 
were dcterniined and presented through computer simulations of the close-loop sys- 
tem. As is evident from the plots, all the suggested control laws were found to take 
the spacecraft to the desired equilibrium point or target attitude through an Eigen- 
axis rotation. The responses of the sustem were found to be of a varying nature 
ranging from highly oscillatory to highly smooth. However only a small section of 
all these possible responses could be presented. A more deeper understanding and 
improvement of the nature of the responses of these control laws may require further 
theoretical analysis of the system and the proposed laws. 




Chapter 5 


Future Work 


In this thesis, it was found that a number of control laws are possible for achieving an 
Eigen-axis rotation from any attitude to any other target attitude. The generaliza- 
tion achieved opens up possibilities for future advancements for better understanding 
and use of the developed control laws. Three main areas where future work can 
concentrate can be identified as folows. 

1. The conditions on <f) and t/i that should be chosen for a particular shape of the 
q and Q trajectory along the Eigen-axis, like non-oscillatory and smooth. This 
requiro.s a chicpcr study into the relation between ^ and ijj and the solutions of 
th<! clo.s<!-l()oi) Control System equations. 

2. Among the diverse range of <l> and ip that can be used, an optimal choice with 
respect to some criteria remains to be found. An optimal-control analysis of 
,the close-loop system with p and ip as inputs is a promising possibility. 

3. A further investigation into the more general and complex form that arose 
during the development of the control law in section 3.1.3 . This may enable 
enhanced control over the trajectory shape within Eigen-axis rotation. 
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Appendix 


A.l Some facts about rotations 

In a 3-dirnontional space, the number of independent variables required to specify the 
position and orientation of a rigid body is 6, which implies 6 degrees of freedom . For 
only rotation, without any translation, the degree of freedom is 3, the remaining 3 
having been lost in fixing the center or origin of rotations. 

Now, lot X bo the position vector of any arbitrary point with respect to some 
rcforeiuu! fraino. For attitude analysis, any rotation of the body can be viewed from 
two frame-works. In the active frame-work, the body is considered to be rotated, 
with the axes or reference frame remaining fixed. In the passive frame-work, the 
body is considered to be fixed, while a reverse rotation of the axes is considered. In 
either case, the above mentioned point has a new position vector say x' which can be 
expressed as a linear transformation of x, given by x' — Ax where A is 3 x 3 matrix. 
For pure rotations, however, the length of the vector must be preserved. Hence, 

x'^x' = x'^x 

=» {AxY'{Ax) = x'^x 

x^A^Ax = x'^x 

=4> x^[A^A-I)x = 0 

As this must be true for any x, we get the important relation 


A^A = I 


(A.1) 
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This is tlu' orthogonality condition, obeyed by any rotational transformation matrix 
A. It is to 1)0 noted that this gives rise to six Equations involving the nine elements 
of the 3 X 3 matrix A, thus leaving only three of its elements totally independent. 
This is consist.ent with the three degrees of freedom for rotations arrived at earlier. 

Now, we consider a reference- frame fixed in inertial space, calling it the inertial 
reference franie, and another body-fixed reference frame attached to the body, Then 
by the attitude of the body we mean the relative orientation of these two frames. 
The attitude of the body for which the two frames coincide, we will call its reference 
attitude. Tlic attitude can be expressed in terms of, say, the angles between the 
axes of tlu’ t.wo frames, or any other similar relevant quantity. The orthogaonal 
matrix A transforms the position vector {xi,X 2 ,X 3 ) of a point on the body in the 
inertial- releriMi<:(' frame to that in the body-fixed frame, {x\,x' 2 ,x'f). Hence, this is a 
representation of the body’s attitude. This transformation matrix equivalently termed 
here as rotation matrix, contains as its elements the cosines of the angles between the 
axes of the two reference-frames, and is thus also called the Direction-Cosine Matrix. 

k 



Figure A.l; The two reference frames 

This can bo shown as follows . We have 

i = [i ■ i')i' + {i ■ j')j' 

= Oiii' + Pik' + jik' 


where 


Q;j = Cosine of angle between i, & i^ 
= Cosine of angle between i.&j' 
7 , = Cosine of angle between i,&k' 


j = {j ■ i')i' -h (j • j')f + U ' k:')k' 
= a2i' + 02 f + 72 ^^ 


Similarly 
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and 


k = {k-i')i' + {k-f)j' + {k.k')k' 

= ot^i' + ^3 -f 'y^k' 

1 lie above gives the relations between the unit- vectors of the two sets of axes. 

But as the vector itself is fixed, we have 

x[i' 4- X 2 j' + Xg/c' = xii 4- x^j + x^k 

= Xiiaii^ + pif + jik') 

-hX2(cX2i’' 4- 02 j' + 72^0 
4-X3(a3i' 4- 03 f 4- 73 /s') 

Xji' 4- Xgj^ 4" XgA:' = (xiGi -}■ X20;2 4- XsOJa)^' 

+ {Xi0i + X202 + X303)f 
+(a:i 7 i 4 - X272 4 - 3:373)^' 


which implies 


x'l 


.•4 


x'3 


ai 0:2 az 
0i 02 0z 
7i 72 73 


Xi 

X2 

3:3 


Evidently, this is = Ax, with A having cosine-elements ai,0i — 

As an example, consider a rotation of the axes through an angle 9 around the 
x-axis, that is, keeping the 2-axis fixed. The Direction-Cosine Matrix for the above 
rotation can be ol)taincd as below. 

Here noting that k and k' are same, and that i,i',j,j' are all perpendicular to k and 
k\ we have 


73 = 1 and a3,/?3,7i>72 = 0 


Also, 


{i ■ i') = (j • f) = cos 9 OL\ = 02 = cos 9 
and [i • j') = {i' ■ j) = cos9 =4> 0i = a 2 = — sin 9 
Hence the Direction-Cosine Matrix for this case is 

cos 9 - sin 0 0 


Az,(? = 


■ sin 6 cos 9 0 

0 0 1 
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Rotationx about, t h<> other axCvS have similar transformation matrices, and successive 
rotations may he n'prescntcxl by successive multiplications by matrices of these type. 
The resutant product matrices, which are the transformation matrices for the resul- 
tant attitudes, have elements as products of trigonometric quantities, and are thus 
numerically difficult to handle. 

This calls for other attitude measures, and we shall study the quaternionic attitude 
measures, which arc based on the important principle of Euler’s theorem. 


A, 2 Euler’s theorem for rigid bodies 

One of the nuxst important facts in rigid body motion, and the basis of the quater- 
nionic rei)rc.s<'ntaion of attitude used in this thesis, is what is well-known in rigid-body 
mechanics as the Euler’s Theorem for rigid-bodies . The theorem is stated as follows. 

Theorem 3 ( Euler’s theorem ) The most general displacement of a rigid-body 
with one point fixed, is a rotation about some axis. 

For a proof of a theorem, two important results should be considered. 

H.esult 1 Any matrix A representing a rigid-body displacement with one point fixed, 
has at least one of its Eigen-values equal to 1 with a corresponding non-trivial Eigen- 
vector. 

Proof: 

To analyse the motion of rigid-body with one point fixed, we consider two reference 
frames, as in previous section. One is is fixed to the body, another being fixed in 
inertial space. In addition, we consider the fixed point on the body as the common 
origin. As before, the attitude may be expressed as a linear transformation between 
the position vectors of a point on the body in the two frames, represented by a; — Ax. 
As the body is rigid, all position vectors have their length preserved. Through an 
analysis exactly as in previous section, A can be shown to be orthogonal. Hence, 

A^A = I 


Now this implies. 


detA^-detA = (detA)^ — 1 

detA = ±1 

It can be shown that the case of -1 is valid for a class of transformations called 
reflections . in which a right-handed axes system is transformed to a left-handed one 
But as we considered both body-fixed and inertial frames to be right-handed, and 
also coincident for the reference rotation, this transformation is not possible as long 
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as the body is rigid- Therefore in this case, the determinant is +1. 
Further, for :i gc'ncrnl orthogonal matrix A, we have 


A^A = I 

{A-I)A^ = I-A^ 

=» dct(j4 — I) • det A^ = det(/ — A^) 

=» det(>l — I) = det(J — A^) 

using dot A‘ = del A = 1. 

Therefore, this implies 

dct(A - J) = det(/ - A) (A.2) 

But from (.lu'ory of (Uiterminants, we have for A being 3x3 

det(/-A) = (-l)'Met(I-A) 

Combining this with Eq A.2 we can write 

det(A-J) = -det(A-/) 

det(A — I) = 0 

From this, it can be concluded that the system of linear equations 


{A -l)x==0 
or Ax — X 

has a non-trivial solution. This evidently implies that A has at least one Eigen-value 
equal to 1 and a corresponding non-trivial Eigen-vector x. 

Result 2 All the Eigen-values of any orthogonal matrix are of magnitude 1 

Proof: Let A be an Eigen-vector of A. A may be complex, or purely real. Hence, 
we have, 

Ax = Air 

It also follows from elementary Linear Algebra that A* is a right Eigen-vector of A^. 
That is, 

= Vi’" 
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Thoiclou., nuilt,ipl\ iiig the hist tAvo cejuations, and using Eq A.l W6 get 

x'\A'^A)x = X*X-x'^x 
=» x^x = A*A • x'^x 


From t.lio al)(m\ il, is evident tluit 


A*A = I A p = 1 

where I A I (k’notos the magnitude or modulus of A. This proves the result. 

Now we arc in a position to give a proof of the Euler’s theorem. Prom the above 
two results, it may be conclu<led that among the three Eigen-vectors, one has a value 
of 1, while tlie other two arc in general complex numbers of magnitude 1. Further, if 
these two are not real, but complex, then they have to be conjugates. This is because 
A being real, its c;hiiracteristic equation 

det(A - XI) =0 

have real coefficients, implying that complex roots occur in conjugates. 

Thus , in general the three Eigen-values can be written as 

A, , 2 ,, = 1 , e® . e-" 

The Eigen- vector corresponding to the Eigen- value 1, is clearly a vector that does not 
get affected by the transformation, suggesting that it may be the axis of a rotation 
transformation. 

The physical significance of the quantity 6 can be understood if we take the three 
Eigen- vectors as our reference-frame, taking the axis Eigen-vector ( that with Eigen- 
value 1 ) {IS the 2 -axis. The resultant transformation A' in this reference-frame is 
given by 


A' = U-^AU 


where U is the matrix containing the 3 Eigen-vectors as its columns, 3rd column be- 
ing the axis-Eigen-vector. It can be seen that A' is the same as the Direction-Cosine 
or Transformation matrix A^^g obtained for the case of a rotation about the 2 -axis, 
through an angle 9, in Sec A.l. Hence, it is clear that any 3x3 Real Orthogonal mar 
trix representing a transformation of a rigid-body with one point fixed is equivalent 
to a rotation about an axis, through an angle, say 6. This is the Euler’s theorem for 
rigid-bodies. 

The axis of the rotation is the Eigen-vactor corresponding to the Eigen-value of 
1, while the other two Eigen- values give the angle of rotation 9, A special 

case is 0 = 0, in which case all Eigen- values are 1, and transformation is an identity 
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tiansfonii.U.ion which aKso can be viewed as a rotation through an angle ^ = 0. An 

CtUiv u<i\ to lind I he angle 0 through which a given A rotates the body can be found 
by observing that 

3 

irace[A] = ^ Ai = 1 + 2 cos 0 (A.3) 

t=i 

From which , 0 is cnisily obtained as 

cos 6 = ^{trace[A] — 1) 

JU 

Fienee, given any r(!al orthogonal matrix A, the axis and angle for the corresponding 
rotation can la; readily found. 


A.3 Quaternions and attitude representation 

At any instant, the attitude of a rigid-body may be considered as resulting from 
a motion, or series of motions, that has taken it from its reference attitude to the 
present one. Reference attitude, it may be recalled, is the attitude for which the body- 
fixed frame coincides with the inertial frame. Now, Euler’s theorem implies that this 
motion or s<!rie.s of motions, is always equivalent to a single rotation about an axis 
through an angle, 'riiis axis, which is fixed with respect to the inertial frame for a 
given altitude, is called the Eigen-axis of that attitude. Therefore the Eigen-axis and 
the angle of rotation can serve to uniquely describe the rotation from the reference 
orientation, tluu'cby describing the present orientation too. Denoting an unit vector 
along the Eigon-axis as n, and the angle of rotation as 9, an algebraic quadruplet 
g = ((? 0 )< 7 i,</ 2 i g:\) called quaternion is defined as below . 

Qq = cos 0/2 (A.4) 

q = h sin 0/2 (A-.S) 

where q = {q\, q-z, (7:i) is considered as the vector part of q, and qo the scalar part. A 
quaternion with scalar part zero is termed a vector quaternion. A norm N{q) of q is 
also defined as 


^{q) = 9o+ 1 q \^— 9o + 9i + ?2 + ql 


It can be checked that for the above definition of q in Eq A.4 and Eq A.5, the norm 
N{g) is always 1. Thus they are precisely called unit quaternions. 

A detailed discussion on the mathematical structure and algebra of quaternions in 
genaral is given in the section A.4. 
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A. 3.1 Relation between rotation matrix and quaternions 

I‘or quatt’nu()us t.o he useful parameters for attitude representation, rotations also 
must he r(q)ros<'ntabl(: by them. This further implies that the attitude transformation 
matrix .4 mu.st be ('xprcssiblc totally in terms of quaternion parameters. Here, the 
relation betweem A and the parameters of q are derived. However, in this case, for 
convcni('nc(\ we chose the active framework. In this frame-work, — it may be recalled, 
the body is ttonsidered to have moved, with the frame fixed. 



Figure A.2: Rotation of vector r about OA 

In the figure A.2, a point C specified by the position vector OC, is rotated about 
an arbitrary axis OA, to the new location C", with position vector OC . The vectors 
OC and OC are denoted by f and f\ while n is a unit vector along the axis OA. 
The angle through which OC rotates is d. 

From the figure, | OA |= (f ■ n), which gives 


II 

O 

(A.6) 

Again 


1 AC 1= f X n 

(A.7) 

Also, from the fig and Eq A.6, 


AO = 00 - OA = f - (n • f) n 

(A.8) 
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TliPioloH’. lioin hqs A. 7 , and A. 8 , and noting that AC and AC are of the same 
niagnilmlt'. wn nan write 


I AC 1— I AC 1=1 f X n 1=1 r — (f ■ h) h\ 

Now, noting the magnitude and direction of AB from the figure, we can write 

AC 


• \ / 


AC 


AB =1 AC I cos^ 

which afi.or using Ec{ A.8 and A.9 becomes 

AC 

AB =1 AC I cos 0 • |- — j = [f - (f • n) n] cos 6 (A.IO) 

Similarly, noting thn magnitude and direction of BC , and using Eq A.9, we obtain 


BC =! AC I sin 6 


f X n 
f X n 


= {f X h) sin 9 


(A.11) 


Finally, wc express the rotated vector f' in terms of f, the axis, and the angle of 
rotation by using the Eqs A. 6, A.IO, and A. 11 as follows 


f = OA + AB + BC 


(f • n) n + [f — (f • n) n] cos 0 + ( 7“ x n ) sin 9 


which after rearrangement gives 

f' = f cos 0 + {f • n)7r[l - cos 0] + ( f x n ) sin 0 (A. 12) 

This is the well-known rotation-formula for finite rotations. At this stage, we intro- 
duce the quaternion, as defined in Eq A.4 and A. 5. Rewriting the above two equations 
for convenience, we have 


qo = cos 9/2 
q = n sin I 

With these definitions, we can expi'ess the terms of the right-hand-side of the rotation 
formula Eq A.12 in terms of the quaternion parameters as below. We have 

r-’cos 0 = f {2 cos^ 0/2 - 1) = f (2qg - 1) = f {ql - q\ - ql - ql) 

{f • n) h [1 - cos 9] = 2 jyf [2 sin^ 9/2] = 2 (5 • f) q 

^ ^ sm y/2 


Similarly, 
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and 


(f X sin 0 = . 2 sin 9/2 • cos 9/2 = 2 (f x ?) ,„ 

With (.lu',s<' suhstit.utions, the rotation formula Eq A. 12 becomes 
r' - r {(/^ -q'^ - q^'j + 2 (q ■ f) g + 2 (f x q) go 


(A.13) 


The above is a relation giving the rotated vector f' in terms of the old vector f and 
the quaternions corresponding to the particular rotation. If the position vectors are 
given in terms of (x,i/,z) co-ordinates with respect to some inertial reference frame, 
then the above relation Eq A.13 can be written as a linear transformation of the 
co-ordinates. This is efvaily obtained by expanding Eq A.13. 



r 



V* 

= 2 








'if + Qo - 5 


— (73?0 

,.2 


QiQa + ? 29 o 
<ii t}2 + gm (ji + 9o - 5 '?2(?3 - Qigo 
'Ml - Mq + ?igo ql + - i 


Xi 

X2 

^3 


(A.14) 


where {x,y,z) and {x',y\z') are the co-ordinates of f and f' respectively, in the 
chosen inertial reference-frame. The above is clearly of the form x' — A{q)x, with 
A{g) being tlie rotation-matrix in terms of g elements. 


A. 3. 2 The reference attitude in quaternions 

The refcreiKX!- attitude has been defined as the one in which the body-fixed frame 
and the inertial frame coincide, and hence, x and x' are same. For this, , obviously, 
A = ly where I is tin; 3x3 identity matrix. Thus, A = I represents the reference 
attitude. 

Let ge be the quaternion representation of the reference attitude, then for 
this, A(g<.) = I. This is satisfied only for the case 

go = 1; gi = g 2 = gs = 0 

Hence, the reference attitude in represented as (1,0, 0,0) in quaternions. This is 
con.sistcnt physically, since this corresponds to the case of 0 = 0 in the defining 
equation.s of g, Eq A. 4 and A. 5, and any n. This just means that the attitude is such 
that the body lias rotated through zero angle about any axes, to arrive at it from the 
reference attitude. A source of ambiguity may be the fact that (-1, 0, 0, 0) also makes 
A(g) = I and hence is also a candidate for the reference attitude representation. This 
reflects the physical fact that the reference attitude may be considered as resulting 
from a 360° rotation from itself, about any axis. This makes $ = 360°, with any n, 
and thus go becomes cos 360° = —1, with g = 0, giving (—1,0,0, 0) as the related 
quaternion. 
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A. 4 Basic quaternion algebra 

Qu;it.(!niiuii.s cousUl.iiU! an algebra by means of which spacial kinematics, especially 
spherical kinematics can be treated elegantly. They were suggested by Hamilton, 
although in a sliglitly different form, but subsequently pushed to the background at 
the advent of vcct.or-algebra. Towards the latter half of this century, interest in it 
was revivc'd lor its suitabilit.y in rotational kinematics. It is closely related to Euler 
Parameter.s, often termed equivalently, and is also similar to spinors in quantum 
mechanics. 

A quat ernion can be defined }U 5 complex number with four units, written as 

Q = qo + qii + q^j + q^k 

where • . q.\ nre real numbers, q^) is called the scalar part and {qi,q2, q-i) the vector 
part. In fact, a ciuaternion with qo = 0 can be considered a 3 -Dimentional Euclidean 
vector, and is c;allcd a vector quaternion. Accordingly, a quaternion can be w'ritten 
formally !us f.ht; sum of a real scalar and a vector. 

q = qo + q 

where q = 91 + 92 + 33 
The addition of quaternions is defined by 

9 + 9' = (90 + qii + 92i + Qzk) + (q'o + 9 i^ + 92i + 93^) 

= (90 + q'o) + (91 + 9l)^ + (92 + 92).? + (93 + 93)^ 

Multiplication of quaternions is based on the following rules for the multiplication of 
the units . 

l.'i = il = i ; l.j = j.l - j \ l-k — k.l = k 

i 2 = f = k'^ = -l 

i.j = —j.i = k ; j.k = —k.j = i ; k.i = —i.k = j 
In addition, it is distributive with respect to summation. With this, we have 

qq' = (90 + 9i^ + Qzj + 93^") (9o + 9l* + 92i + 93^) 

= (9o9o ~ 9i9i ~ 9292 “ 9393 ~ 9494 ) + (9o9i + 9i9o + 9293 ~ 929^^ 

+(9o 92 + 929n + 939i - 9i93)i + (9o93 + 939o + 9i92 “ 929i)^ 
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It, ( ail 1)(' .s('t;u (.hat. inult.iplicat.iou is not commutative. 
post-tiiu!l ipliciUioii can be expressed as a linear operation. 


The premultiplication or 




9i) 

“9i 

“92 

“9.3 

9i 

9o 

-93 

92 

92 

9.3 

90 

~9i 

93 

“92 

9i 

9o 



I 

0 

1 



9l 



to 



CO 

i 



% ~Q2 “93 

9o +93 “92 

92 -93 9o 9'i 

93 92 “9l 9o 



go 


9i 


92 


- ^3 . 


(A.15) 


A cnnjufj<Lt<‘ of a <iuatcrniou (70. 9i, 92, 93) is defined as (go, -gi, -g2, -93). In other 
words, conjugate of 


is given by 


. It follows tliat 


9 = (90, 9) 


9* = (90, “9) 


99* = 9o + 9i + 92 + 93 


This noii-iiogative number is called the ‘Norm’ of the quaternion q, denoted here by 
N{q). Evidently, N{q) = N(g*). 

Quaternioums with N{q) — 1 are called xinit quaternions. 

A. 4.1 Sonic important results for quaternions 

1. Multiplication is associative . q{q'q") = {qq')q" = gg'g" 

2. ( 99 O* = q" ■ 9* 

3. N{qq') = N{q) ■ N{q') 

4. g-i = gVyV(g) 

5. qq' = qtiq'Q - g • g' 4- Qoq' + q'oqq x q' 


Special cases of results 4. and 5. can be written for vector quaternions, as 
4a. g“^ = g* 

5a. qq' = —q • g' + g x g' 
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A. 4. 2 Quaternion multiplications and rotations 

Let be the position vector of a point on the body in inertial frame. . 

For kinematica! analysis, this can be generalized as a vector quaternion X given by 
X = (0,.ti,a:2,.r3). With this, we have the following theorem. 

Theorem A For a vector quaternion X = (0, x) and a unit quaternion q, the opera- 
tion Q\ on X defined by QX = qXq* represents a rotation of x to another vector x' 
in 3- Dim space so that QX = (0, x') 

Proof: Since A" is a vector quaternion, we have A' 4- .AT* = 0 Again, denoting the 
vectorPA' tis A'', wc have X' = qXq*. 

This implies, from the result 2 in subsection A.4.1, 

A" = qX*q* 

Thus, adding the.se two, wc get 

A' + A" =?(A + A*) q* = 0 

Hence, it is clear that X' is also a vector quaternion, equivalent to a 3-dimentional 
vector, expressible as 


A' = QA = (0,x') 

Further, noting that q is unit-quaternion and using result 3, we have 

N(A") = N{ qXq* ) = N{q) • iV(A) ■ N{q*) = N(X) 

which implies x x' — x^ x Hence the above transformation preserves the length of 
the position vector, implying it is a rotation. In fact , by actual multiplication, it 
can be verified that x' and x are related by the relation x' = A{q)x as in Eq A.14. 
Further, if there are two successive rotations, from x to x', and from there to x", 
represented as quaternion operations on A and X' by Q and Q' so that 

X' = QX ; X" = Q'X' = Q'QX 

then we have X' = QX = qXq* 

and A" = Q'X' = g'X'q^ = q'qXq*q' = q"Xq‘ 

where g" = q'q • 

hence, the result Q" = QQ' of two successive rotations is represented by the quater- 
nion q" = qq' which is the product of the quaternions representing the two rotations. 
This is a vei-y useful property, and quaternion multiplication being purely algebraic 
operations, makes numerical processing of rigid-body kinematics much more easier 
and faster. 
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